Berry phase of non-ideal Dirac fermions in topological insulators 
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A distinguishing feature of Dirac fermions is the Berry phase of n associated with their cyclotron 
motions. Since this Berry phase can be experimentaUy assessed by analyzing the Landau-level 
fan diagram of the Shubnikov-de Haas (SdH) oscillations, such an analysis is widely employed in 
recent transport studies of topological insulators to elucidate the Dirac nature of the surface states. 
However, the reported results have usually been unconvincing. Here we show a general scheme for 
describing the phase factor of the SdH oscillations in realistic surface states of topological insulators, 
and demonstrate how one could elucidate the Dirac nature in the real experimental data. 

PACS numbers: 73.25.+i, 73.20.At, 71.70.Di, 72.20.My 



I. INTRODUCTION 



During the last three decades, the Berry phased 
has become an important concept in condensed mat- 
ter physics,^ playing a fundamental role in various phe- 
nomena such as electric polarization, orbital magnetism, 
anomalous Hall effects, etc.- The Berry phase (or ge- 
ometrical phase) in solids is determined by topological 
characteristics of the energy bands in the Brillouin zone 
(BZ) and represents a fundamental property of the sys- 
tem. For example, a non-zero Berry phase, which can be 
measured directly in the magnetotransport experiments, 
reflects the existence of a singularity in the energy bands 
such as a band-contact line in three-dimensional (3D) 
bulk states or a Dirac point in a two-dimensional (2D) 
surface state. ^ Also, the Berry phase of tt is responsi- 
ble for the peculiar "anti-localization" effects in carbon 
nanotubes or graphenei^ Recently, the tt Berry phase has 
been observed in the Shubnikov-de Haas (SdH) oscilla- 
tions in graphenejSiZ' giving one of the key evidences for 
the Dirac nature of quasiparticles in the 2D carbon sheet. 

The 3D topological insulator (TI) also supports spin 
polarized 2D Dirac fermions on its surface?^ which can 
be distinguished from ordinary charge carriers by a non- 
zero Berry phase. Recently, several groups have reported 
observations of the SdH oscillations coming from the 
2D surface states of TIsi^"— In those studies, a finite 
Berry phase has been reported, but it usually deviates 
from the exact tt value. For example, in the new TI 
material Bi2Te2Se (BTS)r^ where a large contribution 
of the surface transport to the total conductivity has 
been observed, the apparent Berry phase extracted from 
the SdH-oscillation data was 0.447r. So far, the Zee- 
man coupling of the spin to the magnetic field has been 
considered''! as a possible source of such a discrepancy. 
Here, we show that in addition to the Zeeman term, the 
deviation of the dispersion relation E(k) from an ideal 
linear dispersion^^ can shift the Berry phase from tt. Wc 
further show how the real experimental data for non- 
ideal Dirac fermions could be understood by taking into 
account those additional factors. 



m/m„=0.15 Vp= 3.4x10 m/s m/mo= 0.125 
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FIG. 1: (Color online) Experimental band dispersions (sym- 
bols) in Bi2Te2Se measured by ARPES in Ref. [H and the 
fitting of Eq. (1) to the surface state (solid line). Large sym- 
bols depict the bulk state. 



II. ENERGY DISPERSION OF SURFACE 
STATES 

The energy dispersion of the surface states in TIs 
can be directly measured in angle-resolved photoemission 
spectroscopy (ARPES) experiments. As an example, Fig. 
1 shows the dispersion of the surface state (together with 
the bulk state) in BTS reported by Xu et al}^ One can 
easily recognize that E{k) is not an ideal Dirac-like dis- 
persion, but it can be fitted reasonably well for the two 
high-symmetry axes with 



E{k) = vphk + /c^ 

2m 



(1) 



with a single Fermi velocity vp — 3.4x10^ m/s and the 
effective mass m which slightly varies with the direction 
in the surface BZ as shown by the solid lines in Fig. 1 
[m/mo = 0.15 (0.125) for the f ^ M (f ^ K) direction 
with mo the free electron mass]. 

Similar fittings can be obtained for other TIs ow- 
ing to the progress in the ARPES studies of these 
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III. BERRY PHASE IN QUANTUM 
OSCILLATIONS 

It is commonly accepted that quantum oscillations 
observed in 3D metals can be well understood within 
Lifshits-Kosevicb^i (the de Haas- van Alphen effect) and 
Adams-Holstein22. (the SdH effect) theories. Recently 
this approach has been generalized to describe magnetic 
oscillations in graphene, which is a 2D system with a 
Dirac-like spectrum of charge carriers i^'^i^'^ There are two 
most prominent features that distinguish such systems 
from materials with a parabolic spectrum: First, rather 
weak magnetic fields are sufficient to bring the system 
into a regime where only a few Landau levels are oc- 
cupied. Second, Dirac quasiparticles acquire the Berry 
phase of tt in the cyclotron motion, changing the phase 
of quantum oscillations. 

In the SdH effect, the oscillating part of p^x follows 

Ap,, ~cos[27r(^-7)], (2) 

where F is the oscillation frequency and 2tt^ is the phase 
factor (0 < 7 < 1). This is the same 7 as in the Onsager's 
semiclassical quantization condition— 

')ttp 

An^—B{N + j), (3) 

when the iV-th Landau level (LL) is crossing the Fermi 
energy Ep {A^ is the area of an electron orbit in the k- 
space) . 7 is directly related to the Berry phase through'' 

,-l = -l-£ndk, (4) 

where il{k)—i J dk u'^i'r) '^k'^^i'^ is the Berry connection, 
M^(r) is the amplitude of the Bloch wave function, F is a 
closed electron orbit (the intersection of the Fermi surface 
E{k) = Ep with the plane kz = const). For spinless 
quasiparticles, it is known^i^^ that the Berry phase is 
zero for a parabolic energy dispersion (7=5) and tt for 
a linear energy dispersion (7 = 0). 

Experimentally, 7 can be obtained from an analysis 
of the Landau-level (LL) fan diagram. There are three 
quantities which are often used as abscissa for plotting a 
LL fan diagram: (i) Landau level index iV, which deter- 
mines the energy E^ of the iV-th LL. (ii) Filling factor 
V {= ^jf^, where Ns is the density of charge carriers, S 

is the area of the sample, = ^ is the number of flux 

quanta, and ^0—^ is the flux quantum), (iii) An integer 
number n which marks the n-th minimum of the oscil- 
lations in px^. Although all three quantities are related 
to each other, the most straightforward way to plot a LL 
fan diagram from the p^x oscillations in a 2D system^ is 



to assign an integer n to a minimum of pxx (or a half- 
integer to a maximum of Pxx)- From Eq. (2), one can see 
that the first minimum in pxx is always in the range of 
< < 1. Thus, the plot of F/Bn vs n, which makes 
a straight line with a unit slope for periodic oscillations, 
is uniquely defined and cuts the n-axis between and 1 
depending on the phase of the oscillations, 7. 

The ordinate 1 / i?„ in a LL fan diagram is determined 
by the Landau quantization of the cyclotron motion of 
electrons in a magnetic field. In 2D systems, upon sweep- 
ing B, Pxx shows a maximum (or a sharp peak in the 
quantum Hall effect^) each time when En{B) crosses the 
Fermi level. Thus, the position of the maximum in pxx 
that corresponds to the A''-th LL, 1/Bn, is given by 

2tt(— -7) = 2nN. (5) 
\Bn I 

On the other hand, the n-th minimum in pxx occurs at 
XjBn when 2n{^ — 7) = 27rn — tt, so the positions of 
the maxima and minima are shifted by ^ on the n-axis. 

The Onsager's relation^ gives F in terms of the Fermi 
wave vector kp as F — {?i/2iTe)TTkp, and this kp can be 
calculated from Eq. (1) as 




Also, when Ep is at the iV-th LL, there is a relation 

ENiBN) = Ep. (7) 
From Eqs. (5)-(7), one obtains 

nu!c \ niVp y mvp I 

where LOc—eB/m is the cyclotron frequency. 

In general case, 7 is a function of -B, meaning that 
oscillations in pxx are quasi-periodic in 1/B. In order to 
calculate 7 one needs to find the eigenvalues Em for a 
given Hamiltonian. 

IV. MODEL HAMILTONIAN 

For the (111) surface state of the Bi2 Sea-family TI 
compounds, the Hamiltonian for non-ideal Dirac quasi- 
particles in perpendicular magnetic fields can be written 

TP 1 

H = Vp{I{xCFy - UyGx) + ^ ^ 2 9s^J■BB<7z, (9) 

where the Landau gauge A = {0,By,0) for the vector 
potential is used, n=?ik+eA, ct^ are the Pauli matrices, 
/is is the Bohr magneton, and gs is the surface g-factor. 
The LL energies are given by^&i^l 

E^^^ = hiu.N ± y 2h vleBN + - 1 g^peB), 

(10) 
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FIG. 2: (Color online) (a) Landau level fan diagram calcu- 
lated for F = 60 T, VF = 3x10^ m/s, = 0, and differ- 
ent m/mo. Arrows show the direction of decreasing m/mo. 
The dashed and dotted lines are the expected behaviors for 
an ideal Dirac dispersion and a parabolic dispersion, respec- 
tively, (b) Landau level fan diagram calculated for F = 60 
T, m/mo = 0.1, ga = 0, and different vf- Arrows show the 
direction of decreasing vf- Insets show the calculated "/{N). 



where and "— " branches are for electrons and holes, 
respectively. The obtained eigenvalues En define the ex- 
act positions of maxima in p^x and, thus, the phase of 
oscillations through Eq. (8). 

In two extreme cases, for non-magnetic fermions (g^ 
= 0), Eq. (8) gives the expected results. First, for a 
linear dispersion (ideal Dirac fermions), m — > cxd leads to 

EN^±^2hevlBN and 7 ^ 5^-^ - TV, giving 7 = 

(Berry phase is it). Second, for a parabolic dispersion, 
VF ^ leads to EN^fiLOc{N + i) and 7 ^ £^ - TV, 

giving 7 = ^ (Berry phase is zero). This gives confidence 
that the expression for 7 given in Eq. (8) is generally 
valid for the topological surface state with a non-ideal 
Dirac cone described by Eq. (1). 
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FIG. 3: (Color online) Landau level fan diagram calculated 
for _F = 60 T and different gs, keeping vf = 3x10''' m/s and 
m/mo = 0.1 constant. Arrows show the direction of changing 
gs- The dotted line is the expected behavior for a parabolic 
dispersion. Inset shows the calculated ^{N). 



V. LANDAU-LEVEL FAN DIAGRAM FOR 
NON-IDEAL DIRAC FERMIONS 

Let us first consider how the LL fan diagram will 
be modified, when both linear and parabolic terms are 
present in the Hamiltonian [Eq. (9)]. For the moment, 
the Zeeman coupHng of the electron spin to the magnetic 
field is assumed to be negligible {qs = 0). Figure 2 (a) 
shows the calculated positions of maxima and minima in 
Pxx for oscillations with = 60 T and vp = 3x10^ m/s 
as m/m^ is varied. One can see that upon decreasing 
m/mQ, the calculated lines on the LL fan diagram are 
gradually shifting upward from the ideal Dirac line that 
crosses the n-axis at exactly ^. Moreover, the lines are 
not straight anymore, which is clearly inferred in the de- 
pendence of 7 vs shown in the inset. With decreasing 
N (increasing B), 7 becomes larger, refiecting the change 
in the phase of oscillations at high fields. 

Similar change in the LL fan diagram occurs if we mod- 
ify another parameter, vp- As shown in Fig. 2 (b), the 
calculated lines are gradually shifting upward from the 
ideal Dirac line as vp is decreased. The results shown 
in Figs. 2 can be understood as a competition between 
linear and quadratic terms in the Hamiltonian [Eq. (9)]. 
Note that for the whole range of the parameters vp and 
m/mo, the positions of maxima and minima in p^x lie 
between two straight lines (shown as dotted and dashed 
lines in Figs. 2) corresponding to 7 = and 7 = ^■ 

Let us now take the Zeeman term into considerations. 
Figure 3 shows the LL fan diagram calculated with F — 
60 T, vp — 3x10^ m/Sj ^-^d m/mo = 0.1, while gs is 
varied. To understand the effect of the Zeeman coupling, 
it is important to recognize the following two points: (i) 
The Zeeman term in Eq. (10) would tend to cancel the 
^hojc term when gs is positive. In fact, when ^TllUc = 
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FIG. 4: (Color online) Landau level fan diagram for oscilla- 
tions in dpxx/dB measured at T = 1.6 K and 9 ~ 0° reported 
in Ref. for BTS. Minima and maxima in dpxx/dB corre- 
spond to n -|- i and n -I- |, respectively. Solid (dark gray) 
line is the calculated diagram for an ideal Dirac cone with vp 
= 3.4x10^ m/s and F = 62 T; dashed (blue) line includes 
the effect of the actual dispersion with m/mo = 0.13; dot- 
ted (red) line further includes the Zeeman effect, where Qs 
= 76 or —45 was determined from a least-square fitting to 
the data. Inset shows the experimental data and calculations 
after subtracting the contribution from an ideal Dirac cone, 

(l/B)Dirac, where A(l/B) = (1/B) - (l/B)Dirac. 

^gsHsB {i.e., gs = 2mo/m) is satisfied, the effect of the 
finite effective mass is canceled and the LL fan diagram 
becomes identical to that for the hnear dispersion (ideal 
Dirac) case. In the present simulations, we use m/mo ~ 
0.1, so that this cancehations occurs when gs — 20. (ii) A 
pair of gs values that give the same l^^iWc — ■^gsf^sBl are 
effectively the same in determining the behavior of the 
LL fan diagram. The resuh of our calculations shown 
in Fig. 3 is a demonstration of these two points. Since 
the Zeeman effect is more pronounced at higher fields, 
the LL fan diagram in Fig. 3 is strongly modified from a 
straight line when the quantum limit is approached, i.e., 
close to TV = 0. 



VI. THE CASE OF BTS 

Let us examine the real data measured in the BTS 
sampler^ in the light of the above considerations. Figure 
4 shows the LL fan diagram for oscillations in dpxx/dB 
measured at T = 1.6 K in niagnetic fields perpendicular 
to the (111) plane In Ref. [ij, the data were simply fit- 
ted with a straight line, and the least-square fitting gave 
a slope of F = 64 T with the intersection of the n-axis at 
0.22±0.12; this result implies a finite Berry phase, but it 
was not exactly equal to tt, which remained a puzzle 
Now, we analyze this LL fan diagram by considering the 
non-ideal Dirac dispersion as well as the Zeeman effect. 
The ARPES data^ for the surface state of BTS (Fig. 



1) gives vp = 3.4x10^ m/s and the averaged effective 
mass m/mo = 0.13. We fix the oscillation frequency F 
at 62 T obtained from the Fourier-transform analysis of 
the dpxx/dB oscillations. "'^^ 

In Fig. 4, the calculated diagram for an ideal Dirac 
cone is shown by the solid (dark gray) line, whereas that 
for the non-ideal Dirac cone with the effective-mass term 
is shown by the dashed (blue) line. One can see that 
the difference is small, which indicates that the effective 
mass of 0.13too is not light enough to significantly al- 
ter the LL fan diagram. One may also see that these 
two lines undershoot the actual data points at smaller n, 
which is even more clearly seen in the inset, where the 
experimental data and the calculations are shown after 
subtracting the contribution from an ideal Dirac cone. 
By further including the Zeeman effect, we can greatly 
improve the analysis, as shown by the dotted (red) line; 
here, gs is taken as the only fitting parameter and a least- 
square fitting to the data was performed. The best value 
of gs is 76 or —45. 

The inset of Fig. 4 makes it clear that it is the slight 
deviation of the experimental points from the ideal Dirac 
line that causes a simple straight-line fitting of the LL 
fan diagram to intersect the n-axis not exactly at 0.5. 
Since the Berry phase in real situations is not a fixed 
value but is dependent on the magnetic field, the simple 
straight-line analysis of the LL fan diagram should not 
be employed for the determination of the Berry phase. 
Obviously, the SdH oscillations of the topological surface 
states are best understood by the analysis which consid- 
ers both the the deviation of the energy spectrum of the 
Dirac-like charge carriers from the ideal linear dispersion 
and their strong coupling with an external magnetic field. 



VII. OTHER MATERIALS 

Similar analysis can be performed for other TIs in 
which the quantum oscillations coming from the 2D topo- 
logical surface states have been observed. Figure 5 shows 
the LL fan diagrams for the SdH oscillations published to 
date for Tl materials^iii'i^ii^'2^ together with the data 
obtained in graphene^ which provides a good reference 
for studies of Dirac fermions. We digitized the pub- 
hshed experimental data in the literature and determined 
ourselves the positions of minima 1/Bmin and maxima 
^/Bmax of the oscillating parts of resistivity (resistance). 
Hall resistivity, or their derivatives with respect to B. 



Material vf (m/s) m/mo Ref. remark 

BiaSes 3.0 xlO^ 0.25 [19| averaged 

BiaTeaSe 3.4 xlO^ 0.13 [12] averaged 

BiaTea 3.7x10^ 3.8 [20] near Dirac point 

graphene 1 x 10® oo [6] calculations 



TABLE I: Parameters of the surface states from ARPES. 
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O BizTej, {dpy,/dB)^i„ 

• BijSej, (AR,,)„i„ 

o graphene, (R„)^in „,3^ 

■ Bi2Te3 nanoribbon, 

(^''xx)min, max 
□ BTS, (dp,^/dB)^i„ ,^3^ 



FIG. 5: (Color online) Landau level fan diagrams for SdH 
oscillations observed in various TIs and graphene. Symbols 
are obtained from the published experimental data in the 
literature. Solid lines are calculations taking into account 
the non-ideal dispersions of the surface states (determined 
by m/mo) and the Zeeman coupling to an external magnetic 
field (determined by Qs)- Dashed lines are calculations for 
ideal Dirac fermions (m/mo = oo and Qs =0). Open dia- 
monds are {dpy^ / dB)min,max in Bi2Te3 from Ref. H; filled 
circles are {/!^Rxx)min in Bi2Se3 from Ref. [Til : open circles 
are {Rxx)min,max in graphene from Ref. filled squares are 
{/^Rxx)min,m,ax iu a Bi2Te3 uauoribbon from R ef. Ila : open 
squares are {dpxx /dB)min.max in BTS from Ref. [12. 



Material 


Ref F (T) Ef (cV) 




9s 




BiaSeg 


[in 


30.7 


0.074 


55 


or 


-39 


BiaSeg 




88.6 


0.143 


55 


or 


-39 


BiaTeaSe 


[12] 


62.0 


0.152 


76 


or 


-45 


BiaTeg 


[9] 


27.3 


0.074 


65 


or 


-65 


BiaTea, nanoribbon 


[16] 


54.7 


0.101 


65 


or 


-65 


graphene 


m 


43.3 


0.239 










TABLE IL Parameters used for the calculations shown in Fig 
5. 



The obtained data for various materials are plotted as 
functions of n in Fig. 5. Note that, to avoid ambiguities, 
we considered only those data that show oscillations with 
a single frequency<^ 

The parameters of the surface states used in our fan- 
diagram analyses have been obtained from the published 
ARPES data by fitting them in the same way as for 
BTS (see Fig. 1). Tabic I shows vp and mjmQ for the 
Bi2Se3/Bi2Te3 family and graphene. These parameters 



were fixed during the fitting of the data shown in Fig. 5. 
The only parameter that could vary in our calculations 
was Qs- Note that the frequency of oscillations F (and, 
thus, the Fermi energy Ep) is essentially determined by 
the periodicity of the observed oscillations. Table II sum- 
marizes the parameters thus obtained. The results of our 
calculations are shown in Fig. 5 by solid lines. Dashed 
lines depict the behavior expected for ideal Dirac cones 
(m/mo=oo) and negligible Zeeman coupling (g^. = 0) for 
the TI data . One can clearly see in Fig. 5 that only 
graphene shows the ideal behavior in the LL fan dia- 
gram: a straight line that crosses the n-axis at 0.5. All 
TI materials, despite their essentially Dirac-like nature 
of the surface state, present the LL fan diagrams that 
deviate from the ideal behavior. (The deviations from 
the dashed lines are most clearly seen in strong magnetic 
fields.) 

In view of the good agreements between the data and 
the fittings for all the materials analyzed in Fig. 5, 
one may conclude that the advanced analysis consider- 
ing both the curvature of the Dirac cone and the Zeeman 
effect can reasonably describe the SdH-oscillation data 
obtained for TIs and confirm the Dirac nature in their 
surface states. 



VIII. SUMMARY 

We derived the formula for the phase 7 of the SdH 
oscillations coming from the surface Dirac fermions of 
realistic topological insulators with a non-ideal disper- 
sion given by Eq. (1). We also calculated how the cur- 
vature in the dispersion as well as the effect of Zeeman 
coupling affect the Landau-level fan diagram of the SdH 
oscillations for realistic parameters. Finally, we demon- 
strate that the Landau-level fan diagrams obtained from 
recently reported SdH oscillations in topological insula- 
tors can actually be understood to signify the essentially 
Dirac nature of the surface states, along with a relatively 
large Zeeman effect in those narrow-gap materials. 
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